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Abstract. Using the Anyon-Hubbard Hamiltonian, we analyze the ground-state 
properties of anyons in a one-dimensional lattice. To this end we map the hopping 
dynamics of correlated anyons to an occupation-dependent hopping Bose-Hubbard 
model using the fractional Jordan-Wigner transformation. In particular, we calculate 
the quasi-momentum distribution of anyons, which interpolates between Bose-Einstein 
and Fermi-Dirac statistics. Analytically, we apply a modified Gutzwiller mean-field 
approach, which goes beyond a classical one by including the influence of the fractional 
phase of anyons within the many-body wavefunction. Numerically, we use the density- 
matrix renormalization group by relying on the ansatz of matrix product states. As a 
result it turns out that the anyonic quasi-momentum distribution reveals both a peak- 
shift and an asymmetry which mainly originates from the nonlocal string property. In 
addition, we determine the corresponding quasi-momentum distribution of the Jordan- 
Wigner transformed bosons, where, in contrast to the hard-core case, we also observe 
an asymmetry for the soft-core case, which strongly depends on the particle number 
density. 
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1. Introduction 

A fundamental principle of quantum statistical mechanics in three dimensions is the 
existence of two types of particles: bosons obeying Bose-Einstein statistics and fermions 
obeying Fermi-Dirac statistics. However, in a two-dimensional electron liquid, quasi¬ 
particles made of electrons in the fractional quantum Hall effect are charged anyons 
obeying fractional statistics [1-5]. While exchanging two anyons, the many-body 
wavefunction acquires a fractional phase where the statistical parameter 6 varies 
with the interval 0 < 6 < n and corresponds to a fractional statistics. As compared to 
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bosons and fermions, anyons exhibit a wide range of previously unexpected properties 
and the concept of anyons plays an important role in numerous studies of condensed 
matter physics and of topological quantum computation [6-13]. In order to generalize 
the Bose-Einstein and Fermi-Dirac statistics by allowing a maximal hnite integer particle 
number occupying the same quantum state, Gentile supplied an intermediate statistics 
that had been proven to be also valid for a g-fermion [14-19]. However, Y. Shen et 
al. pointed out that the anyon statistics was not a complete analogue of the Gentile 
statistics [20]. Based on a generalized Pauli exclusion principle, Haldane provided a 
useful concept of fractional statistics in arbitrary dimensions [21]. Polychronakos also 
suggested another form of the fractional exclusion statistics [22]. However, it was shown 
via the virial expansion that all these fractional statistics proposed in the literature do 
not apply for anyons [20,23]. 

In one dimension (ID), anyons were realized as low-energy elementary excitations 
of the Hubbard model of fermions with correlated hopping processes [24]. Alternatively, 
it was suggested to create anyons by bosons with complex-valued occupation-dependent 
hopping amplitudes by photon-assisted tunneling in ID optical lattices [25,26]. Recently, 

S. Greschner and L. Santos proposed a Raman scheme to improve the proposal of 

T. Keilmann et al. in Ref. [25] and deduced a rich ground-state physics including 
Mott-insulators with attractive interactions, pair-superfluids, dimer phases, and multi- 
critical points [27]. An even simpler scheme for realizing the physics of ID anyons with 
ultracold bosons in an optical lattice has recently been proposed in Ref. [28]. It relies 
on lattice-shaking-induced resonant tunneling against the energy off-sets created by the 
combination of both a potential tilt and on-site repulsion. In contrast to the above 
mentioned former proposals based on internal atomic degrees of freedom, no lasers 
additional to those already used for the creation of the optical lattice are required [28]. 

The physical properties of ID anyons are intriguing and complicated. Theoretically, 
using the generalized coordinate Bethe ansatz method, Kundu obtained the exact 
solution of a ID anyon gas with the delta-function potential, which shows that the 
effective interaction among the anyons is modified by the anyonic statistical parameter 
[29]. Furthermore, the ID anyon gas with a special interaction potential has been 
investigated in various ID systems in order to get the ground-state properties of the one¬ 
dimensional anyon gas [30-35]. For anyons in the ID lattice systems one investigated, for 
instance, the statistically induced phase transition [25], quantum scaling properties [24], 
dynamics properties of hard-core anyons [36,37], quantum walks [38,39], and anyonic 
Bloch oscillations [40]. For the hard-core anyons in the ID optical lattice, Y. Hao et 
al. investigated how the fractional statistics affected the ground-state properties by 
mapping the hard-core anyonic Hamiltonian to a noninteracting fermionic system with 
a generalized Jordan-Wigner transformation [41]. By calculating the one-particle Greens 
function of the ground state, the reduced one-body density matrix and, thus, the quasi¬ 
momentum distributions can be obtained for different statistical parameters of anyons. 
The results showed that the momentum distributions in the Bose and Fermi limit turned 
out to be symmetric, but those of anyons were in general asymmetric and shifted due 
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to the fractional statistics that anyons obey. 

In this paper we complement the previous studies of the ground-state properties 
of a ID quantum gas of anyons conhned in optical lattices. To this end we map the 
Anyon-Hubbard model to an occupation-dependent hopping Bose-Hubbard model with 
the help of a fractional version of the Jordan-Wigner transformation. This mapping has 
the consequence that the Hilbert space of anyons can be constructed from that of bosons, 
so that one has access to the two-point correlation function of either the original anyonic 
or the Jordan-Wigner transformed bosonic creation and annihilation operators. With 
this we investigate the quasi-momentum distributions of either bosons or anyons, the 
latter interpolating between Bose-Einstein statistics and Fermi-Dirac statistics. Firstly, 
in the hard-core limit, we determine the quasi-momentum distributions of anyons with 
density-matrix renormalization group calculations, which numerically reproduces the 
results of Y. Hao et al. [41]. In addition, by suggesting a modihed Gutzwiller mean-held 
approach to include the inhuence of the fractional phase of anyons within the many-body 
wavefunction, we obtain an approximative analytic expression for the numerical results 
of Ref. [41] both for a hnite system and in the thermodynamic limit. In particular, we 
analyze in detail how peak-shift and asymmetry of the quasi-momentum distribution 
of anyons depend on both the fractional phase 6 and the particle number density uq. 
Furthermore, we extend the hndings of Ref. [41] by working out also the more general 
soft-core case. The quasi-momentum distribution of anyons reveals the pseudofermion 
property at the fractional phase 6 = n. Surprisingly, the quasi-momentum distributions 
of the Jordan-Wigner transformed bosons shows a density-dependent asymmetry which 
is not found for the hard-core case. 

The outline of this paper is as follows: in Sec. 2 the Anyon-Hubbard model and the 
mapping between anyons and bosons are discussed. Both the classical and the modihed 
Gutzwiller mean-held approach is introduced in Sec. 3. The ground-state properties 
of the ID Anyon-Hubbard model are determined by studying the quasi-momentum 
distribution of anyons and bosons in both the hard-core and the soft-core case in Sec. 
4 and Sec. 5 respectively. Gonclusions are given in Section 6. 

2. Anyon-Hubbard model 

The hopping dynamics of correlated anyons on a ID lattice is described by the Anyon- 
Hubbard (AH) Hamiltonian [25] 

TT ^ 

= -J + h.c.) + - 1), (1) 

i=i 1=1 

where J > 0 denotes the tunneling amplitude connecting two neighbouring sites, U 
stands for the on-site interaction energy, fij = aja,- represents the number operator at 

^ t ^ 

site j, and the operators aj, aj create or annihilate an anyon on site j. For ID anyons. 
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the operators and dj obey the generalized commutation relations [25,29] 

a,.a* - , 

djdk - = 0 , 


where 9 is the statistical exchange phase, the sign function sgn(j — k) = —1,0,+1 for 
j < k, j = k, and j > k, respectively. Note that, since the sign function sgn(j — k) = 0 
for j = k, two particles on the same site behave as ordinary bosons irrespective of the 
statistical parameter 9. Moreover, anyons with the statistical exchange phase 9 = n are 
pseudofermions, i.e. while being bosons on-site, they are fermions off-site. 

Here we use an exact Anyon-Boson mapping in ID in terms of a fractional version 
of the Jordan-Wigner transformation [25] 


dj = bj exp 



( 3 ) 


where the number operator reads fii = ajoj = blbi and hi, hi are bosons operators 
following the commutation relation [5*, 6]] = 6ij and [&|,&]] = 0 = [hi,hj]. Inserting the 
Anyon-Boson mapping (3) into Eq. (1), the Hamiltonian can be rewritten as [25] 




Wh- 


u 




hj{hj - 1 ), 


( 4 ) 


t=i 


where the conditional hopping of bosons from right to left, i.e. j -|-1 —?■ j, occurs with an 
occupation-dependent amplitude If the target site j is unoccupied, the hopping 

amplitude is simply J. If it is occupied by one boson, the amplitude becomes complex 
and reads , and so on. We emphasize that the hard-core limit of anyons in Eq. (1) 
coincides with that of bosons in Eq. (4) due to hi = d\di = h\hi. It is also clear that 
the mapped bosonic Hamiltonian (4) describes local occupation numbers beyond the 
hard-core limit rij > 1. 

Using the nonlocal exact mapping (3) between anyons and bosons, the Anyon- 
Hubbard Hamiltonian (1) leads to a bosonic Hamiltonian Eq. (4) that can be solved 
either analytically or numerically in order to determine the ground-state properties of 
anyons in ID lattice systems. Theoretically, since the reflection parity symmetry in the 
Hamiltonian is broken, we suggest below a modihed Gutzwiller mean-held, which goes 
beyond the classical one found in literature [42-45] in order to include the inhuence 
of the fractional phase of anyons on the hopping dynamics. To get numerical results, 
we use the density-matrix renormalization group (DMRG) [46-50], which was already 
applied to the realm of anyons, for instance, in Ref. [25]. To this end, we rely on the 
ansatz of matrix product states (MBS) with the system length L and open boundary 
conditions, which is more efficient than periodic boundaries [51-55]. The code is based 
on a variational ansatz using MBS with the restricted subspace of integer hlling, where 
our simulations admit a maximum of hve particles per site and the maximum bond 
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dimension of MPS equals to 1000 [56]. Moreover, Refs. [53,54,56] show the respective 
details how to calculate the expectation of the string operator in Eq. (3) and the 
correlation functions in DMRG with MPS ansatz. 


3. Mean-field approximation 


At hrst we work out a Gutzwiller mean-held approach, which turns out to provide 
qualitative satisfactory results for the quasi-momentum distributions. In a Gutzwiller 
(GW) approach for bosons, the many-particle state |G) is generically approximated by a 
product state of single lattice-site states 1$^), which can be expressed as a superposition 
of different number states on a lattice site [45] 

( "^max \ 

(5) 

n=0 J 

Here Umax represents a truncation at some sufficiently large maximal particle number 
and fn'^ denotes the Gutzwiller probability amplitude of hnding n bosons on site j, 
which is normalized such that 

^max 

= ( 6 ) 

n=0 

The total energy of the system 


£tot — 


int 


£-1, + £\ 


kin 


(7) 


then follows from determining the ground-state expectation values of the respective 
operators according to the Appendix. As a consequence, the interaction energy (A.l), 
the chemical potential term (A.2), and the kinetic energy (A.3) depend on the Gutzwiller 
coefficients fn \ Due to the polar decomposition fn^ = Fn^e^°‘’^ , the total energy (7) 
has to be minimized with respect to both the absolute values Fn'^ and the phases 

In the classical GW mean-held approach of a homogeneous system, to which we 
refer as cGW in the following, one assumes that the ground state is a product of identical 
states on each of the L lattice sites [42-45], i.e., 

|t;> = |t>. = (8) 

This has the consequence 

A’ = U. . (9) 

i.e. one has the identical wave function on each lattice site. 

Due to the occupation-dependent amplitude for nearest-neighbour hopping 

in the bosonic Hamiltonian (4), the rehection parity symmetry is broken. Therefore, we 
suggest to drop the cGW constraint of having identical states on each site to a modihed 
GW mean-held approach (mGW). To this end, we note that the kinetic energy Eq. 
(A.5) in the Appendix depends on the diherence /3n ^ = On ^ of the phases of the 
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probability amplitudes fn'^ and of finding n and n +1 bosons on site j, respectively. 
If all phase differences would be identical, we would recover cGW. Therefore, we dehne 
mGW according to 

F«) = F„, A/3P =/3P-/?«+')= A/J„, (10) 

In order to fix the values of the phase differences fin \ we complement (10) by the 
additional assumption 

= 2/3„ . (11) 

This means that the absolute value Fn'^ of the probability amplitude fn'^ is identical on 
each site and the non-vanishing phase difference ^j3n^ between two nearest-neighbour 
sites is the same in the whole chain. 

The application of cGW and mGW is performed in the Appendix for both the 
hard-core case in Sec. Appendix A.l, which means Umax = 1, and the soft-core case by 
assuming Umax = 2 in Sec. Appendix A. 2. Provided that Gutzwiller amplitudes and 
phases have been determined, one can calculate the two-point correlation function of 
bosons, which has the rather simple expression 

(6l6,) = i„(ft.) + (l-iy)(6.')(6,). (12) 

The correlation function of anyons can be written as 

{a\aj) = 5ij{hi) + (1 - . ( 13 ) 

For the case i < j and i > j we get correspondingly. 



(14a) 

i<l<j 


4^,) (11 

(14b) 


j<l<i 


where the respective expectation values depend on the Gutzwiller amplitudes fn'^ 
according to Eq. (A.4) in the Appendix. 

Furthermore, in order to investigate the effect of the anyon statistical parameter 6, 
we will focus on density distributions of bosons and anyons in the quasi-momentum 
space, which are believed to be interesting [25] and are dehned via the Fourier 
transformation of the correlation function 



(15a) 


(15b) 


ij 


In the following, we are going to hnd the ground state of a homogeneous system 
for a given particle number density uq > 0, ratio J/U, and varying statistical parameter 
6. Thus, combining the analytic approaches cGW, mGW with numerical results 
from DMRG, we determine the correlation functions and, thus, investigate the quasi- 
momentum distributions of the ground state. 
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4. Hard-core anyons 

In order to illustrate the applicability of the Gutzwiller mean-held theory, we investigate 
at hrst the simplest case of hard-core anyons. This means that we assume U/ J —)■ cxo, 
where bosons are impenetrable and each site contains at most one particle, i.e., rimax = 1- 
In the following, we present only the calculations of the correlation function and the 
quasi-momentum distribution, while the mathematical details for minimizing the total 
energy are relegated to Sec. Appendix A.l of the Appendix. 

According to Eq. (12), the correlation function of bosons reads 

= dijUQ + (1 - (5p)no(l - , (16) 

which reduces to 

= 5ijnl + uo(l - no), 

since A/Sq = 13^'^ — = 2m7r (see Sec. Appendix A.l). Using Eq 

ground-state quasi-momentum distribution of bosons can be written as 

= n^-Fno(l-no)L(5fc,o, (18) 

which is ^-independent and the number of particles in the condensate is Nq = A(1 — no). 
Fig. 1 shows the resulting ground-state quasi-momentum distribution in a hard-core 
system at density no = 0.5. 

Since in Eq. (18) is 6'-independent in the GW approach, if one changes 6 from 

0 to TT, all quasi-momentum distributions with different 9 coincide as shown in Fig. 1 
(a). The quasi-momentum distribution of bosons in DMRG as shown in Fig. 1(b) 

turns out to be also ^-independent because the phase factor disappears in the bosonic 
Hamiltonian (4) in the hard-core limit and the systems reduces to a hard-core Bose- 
Hubbard model. As a consequence, bosons condense at zero quasi-momentum in both 
the GW approach and DMRG. However, the bosonic population N{1 — hq) of fc = 0 in 
the GW approach is much larger than that in DMRG because the mean-held approach 
does not contain quantum huctuations which generally broaden the quasi-momentum 
distributions. 

In the Bose limit 0 = 0, we have dj = bj and {d\dj) = (g&j), thus 

= 0)) = (fi'"’). (19) 

which is shown in Fig. 1(a) and (b) for the density uq = 0.5. 

According to Eqs. (14) and (A.4d) in the Appendix, the correlation function of 
anyons reduces to 

{d\dj) ^ no(l - no) [(1 - no) + noe'®] , (20a) 

i<l<j 

{d\dj) ^ no(l - no) [(1 - no) + noe"'®] . 

j<l<i 


(17) 

(15a), the 


(20b) 





Figure 1. Quasi-momentum distribution of the ground state in a hard-core system 
(no = N/L = 0.5, L = 120). (a) Gutzwiller approach and (b) DMRG for anyons and 
bosons. 


Using Eq. (15b), the ground-state quasi-momentum distribution of anyons can be 
written in the ^-dependent form for a hnite open system 


L-l 


= no+no{l-no)-\{L-l)(e 




—ikm m—lj^ikm *m—1 
c Z Z 


( 21 ) 


m=2 


where 2 ; = = x + iy = (1 — uq) -f , thus r = 1 if 6^ = 0 or uq = 1, otherwise 

r < 1. 

In the Fermi limit 6 = tt, we obtain from (21) 



no -b 2no(l - uq) 


cos(/c) — X 
1 — 2x cos(/c) -b x‘^ 


(1 -a:-^)[cos(/t)(l + a;^) - 23 ;] ) 
L[1 — 2a;cos(/t)-b /’ 


where (h^) = a: = 1 — 2no G (—1,-bl) if 0 < no < 1. In the 

thermodynamic limit, this reduces to 




{n'J) no -b 2no(l - no) 

|x|<l 


cos(/c) — X 


1 — 23:cos(fc) -b x'^ 


(23) 


For the system at no = 0.5, the GW approach gives 

(n^P) = 0.5 + 0 . 5(1 - y) cos(k) 

1j 

(nb) — 0.5 + 0.5cos(l:), 


( 24 ) 
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which is shown as the black solid line in Fig. 1(a). However, in Fig. 1(b) from DMRG, 
the quasi-momentum distribution of anyons in the Fermi limit represents a step-like 
distribution, which is the characteristic feature of free fermions. This shows that DMRG 
can grasp the fermion-like feature of anyons with statistics 6 = tt, whereas the GW 
approach achieves this only approximately. 

If iFmax^ defines the quasi-momentum, where has its maximum, we have 

Rmax(6* = 0) = 0 and K^IA O = tt) = 0, i.e., both the quasi-momentum distribution 
of the Bose and of the Fermi limit are symmetric about zero quasi-momentum. Indeed 
this results from 5^ 0 in Eqs. (18) and (19) in the Bose limit and from calculating 
d{n^l'’)/dk = 0 of Eqs. (23) and (24) in the Fermi limit. 

In the fractional phase interpolating between the Bose and Fermi limit 0 < 6^ < vr, 
we find from (21) 



no(l — rio) r(T —1) —L(x;e ^^)-\-{ze 


■4-0.c. 


(25) 


In the thermodynamic limit, this reduces to 



no-F2?7,o(l-no) 

r<l 


cos(/c) — r cos(<yc) 

1 — 2r cos(/c —</9)-fr^’ 


(26) 


where r = |1 — no 4- and = arg(l — no -|- noc'®). Fig. 1 clearly show that: i) the 
respective peak of the quasi-momentum distribution is shifted to positive momentum, 
ii) the quasi-momentum distribution is asymmetric about its own peak, and iii) the 
quasi-momentum distributions broaden and flatten with increasing 9. Gomparing Fig. 
1 (a) with (b) or (p^hj) in Eq. (17) with {a\aj) in Eq. (20), we conclude that the 
nonlocal string property of anyons, which is denoted by ~ ’^o) + or 

nj<z<i[(l ~ ’^o) + rioe“^®] in Eq. (20), is the reason why anyons prefer to appear on 
positive momentum. This is explained numerically in Ref. [41]. 

In order to find the maximum iFmax of the quasi-momentum distribution, we 
analytically evaluate d{fp^)/dk = 0 for Eq. (26), which gives 


= 


TT— /C, uq > 2(\/2 —1) and "di < 6'< 792, 
/C, Otherwise, 


(27) 


where we have introduced the abbreviations 


n ^ • /2 — TT-Ox „ StT „ 

Vi = J + arcsm( ), ^^2 = y - W , 

. sin6'(2no — nn)-|-nnsin6*cos6' 

/C = arcsm --—-- -tt---— 

(1 — cos6*)(l — no)^ + (1 + cos^) 


(28) 

(29) 


Fig. 2(a) shows the resulting fractional phase ^-dependence of anyonic iFmlx at 
density hq = 0.75, 0.5, and 0.25. The DMRG results (L = 120) clearly show that the 
peaks of the quasi-momentum distributions increase linearly from Kmlx{0 = 0) = 0 
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Figure 2. (a) Maximum of quasi-momentum plotted against the 

fractional phase 9 of anyons at density uq resulting from DMRG (L = 120) and GW 
in the thermodynamic limit L —>■ oo. (b) it'max plotted against the inverse system size 
at density no = 0.75 resulting from GW. (c) Maximum of the finite system size 

L divided by that of the infinite system size L —>■ oo plotted against the inverse system 
size at density ng = 0.75 resulting from GW. 


(bosons) and come back to = vr) = 0 (fermions) while increasing 6. The GW 

approach shows the same linearly increasing property at small 9, i.e. near bosons. This 
can be quantitatively understood by taking the limit 6* —)■ 0 in Eqs. (27) and (29): 

K ^l^ {9 —)■ 0) ~ arcsin(no6*) ~ tlqO . (30) 


At the Fermi side 6 ^ —)■ tt, we have correspondingly 




~ arcsin 


ripij - 9) 
1 - no 


nojn - 9) 
1 - no 


(31) 


where the slope of i^mix from GW is smaller than that from DMRG. In general, the 
GW approach always yields a smaller shift of positive momentum in comparison with 
the DMRG results for any statistical parameter of anyons between the Fermi and the 
Bose case. Moreover, we hnd that both DMRG and GW give a larger peak shift for a 
bigger particle number density no- 

In order to determine the thermodynamic limit of the quasi-momentum 
distribution. Fig. 2 (b) and (c) show the system size dependence of and 
at density no = 0.75 in the GW approach, where is the maximum value of the quasi¬ 
momentum distribution with the hnite system size L and corresponds to the 
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thermodynamic limit. Considering the uncertainty of the quasi-momentum in the hnite 
system Ak = 27r/L, iCmax turns out to be independent of the system size. According 
to Fig. 1, the bosons condense at zero momentum with a sharp distribution and the 
fermions have a flat distribution. In Fig. 2(c), with 6 increasing from the Bose to the 
Fermi case, we hnd that the condensate population near the Bose side slightly depends 
on the system size, but the maximum particle number near the Fermi side is almost 
independent on the system size, since the behavior of anyons with 0 ~ tt looks much 
more like fermions without condensation. Note that we do not show here the case of 
bosons 6 = 0 because bosons condense at zero momentum. 


5. Soft-core anyons 

Before going to the soft-core case, we remind that within the Gutzwiller theory (b\bj) 
and in Eqs. (17) and (18) turn out to be independent of the fractional phase 

9 of hard-core anyons because A^q = 2m7i and the local state |2) is not occupied, 
i.e. F 2 = 0. Considering the occupation-dependent hopping in the soft-core case, we 
are going to discuss the influence of A/3i and non-vanishing F 2 in the following. The 
respective details of minimizing the total energy are given in Sec. Appendix A.2 of the 
Appendix. 


5.1. Boson Quasi-momentum Distribution 

According to Eq. (12) and supposing Umax = 2, the correlation function of bosons is 

(bib,) = 6,,no + (1 - S,,)[A + , (32) 


where A = F^^(F^ + V 2 F 0 F 2 ), B = Ff{V2FoF2 + 2F|) and A + 5 = C. Using Eq. 
(15a), the quasi-momentum distribution of bosons can be written as 



no — C + ALbkfl 


1 - cos[(fc + g)L] 
L[1 — cos(/c -|- 9)] 


and yields in the thermodynamic limit 


(33) 


{nf'>) no-C + A5{k) + B5{k + 6). (34) 

Here, compared with the 6'-independent of the hard-core case in Eq. (18), the 

quasi-momentum distribution of bosons from mCW has a background distribution uq—C 
and two peaks at A; = 0 and k = —9 with the intensities A and B, respectively. 

Noticeably, if comparing the intensity A = F^(Fq -h '/ 2 F 0 F 2 ) of the k = 0 peak with 
the intensity B = F^{\/2FoF2 -\- 2F|) of the k = —9 peak, we hnd that the discrepancy 
between the intensities A and B is determined by the occupation F^ of |0) and Ff 
of |2). As these occupations Fq and F 2 change with the density uq, the resulting 
quasi-momentum distribution of bosons reveals a characteristic density dependency. 
Therefore, we discuss now different density regimes. 
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quasi-momentum k 

Figure 3. Quasi-momentum distribution of bosons with different statistical phase 0 
{J/U = 0.1, L = 120). (a), (b) and (c) from mGW at density uq = 0.25, uq = 1.25, 
and no = 0.75. (d), (e) and (f) from DMRG at density no = 0.25, no = 1.25, and 
no = 0.75. 

In the dilute limit no 1, the positive Hubbard U depresses the state |2) and 
F 2 ~ 0, thus the properties of are similar to that of the hard-core case. Fig. 3(a) 

shows the quasi-momentum distribution of bosons at density uq = 0.25 and J/U = 0.1 
in the mGW approach. For 6* = 0, we have = no — C + C6{k), which demonstrates 

that there exists only one peak of the intensity C = A + Beitk = 0 and the quasi¬ 
momentum distribution is symmetric. For the case of 0 < 0 < vr, we have (34), which 
demonstrates that there exists one much lower peak of the intensity B at the negative 
momentum k = —6 than that of the intensity H at fc = 0 due to <C Fq. Here, 
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we should note that ~ '^))i which is plotted for slightly asymmetric interval 

—TT < A; < TT in Fig. 3(a), is also a symmetric distribution since has period 271. 

The results of DMRG with the same parameters are shown in Fig. 3(d), where the 
symmetric quasi-momentum distribution of 0 = 0 is a little sharper than that of 6^ = vr. 
But for the case of 0 < 6* < tt, there exists an asymmetric quasi-momentum distribution 
and more bosons have negative momentum according to the inset of Fig. 3(d). In 
general, quantum fluctuations broaden from mGW. Taking this into account, 

if we would add the effect of quantum fluctuations on the mGW quasi-momentum 
distribution, the peak at /c = 0 and k = —9 would broaden. Since the peak at 
k = —9 from mGW is pretty low, a broadening would be merged into the background 
of the quasi-momentum distribution. Therefore, at small particle densities, the mGW 
approach gives qualitatively the same quasi-momentum distribution as DMRG, i.e. the 
peaks locate at A; = 0 and more bosons have negative momentum. 

For the high density uq = 1.25 in Fig. 3(b) and (e), there exists more than one 
particle per site on average. Thus the vacuum state |0) is almost unoccupied in the 
ground state, i.e. Fq <C Ff. For 0 < 9 < tt, Fig. 3(b) clearly shows that the high peaks 
of the quasi-momentum distribution exactly appear aX k = —9 in comparison with that 
at /c = 0. The results of mGW qualitatively coincide with that of DMRG in Fig. 3(d) 
although the momentum value k of the DMRG peaks is a little larger than —9 because 
the peaks aX k = —9 compete with the peaks at A: = 0 under the influence of quantum 
fluctuation. Of course, 6^ = 0 and 9 = tt still give the symmetric quasi-momentum 
distribution of bosonic operators. 

Fig. 3(c) and (f) show the quasi-momentum distribution at the density tiq = 0.75 
in the middle of uq = 0.25 and uq = 1.25. The competition between the A: = 0 peak of 
the intensity A and the k = —9 of the intensity B of the quasi-momentum distribution 
leads to hnal peak shifts from A; = 0 to the negative momentum but comes back to A; = 0 
while increasing 9 from 0 to tt as shown in Fig. 3(f). Note that this peculiar behaviour 
is not reproduced by mGW in Fig. 3(c). 

If we think about the occupation-dependent factor in (4) again, we recognize 
in Fig. 3 a clear physical picture of the ground state: i) for a small density uq < 1, 
the states |0) and |1) are preferred and A/Xq = 2m7r is the leading term, which gives 
a quasi-momentum peak of bosons at A; ~ 0; ii) for a density 1 < uq < 2, the states 
|1) and |2) are preferred and A/Si = 2l7r — 9 is the leading term, which gives a quasi- 
momentum peak of bosons at A; ~ —9. In order to deal with the density 2 < uq, we 
should truncate the mGW approach at Umax = 3,4, • ■ •. For example assuming rimax = 3, 
the quasi-momentum distribution of bosons shows an additional peak of A: = —29. For 
the density 2 < uq < 3, the leading peak of the quasi-momentum distribution of bosons 
appears at the momentum k ~ —29. 
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hard-core 


soft-core 



quasi-momentum k 


Figure 4. Quasi-momentum distribution of anyons with different statistical phase 0 
{J/U = 0.1). (a), (b), (c) and (d) from DMRG {L = 120) at density uq = 0.25 and 
no = 0.75. (a'), (b'), (c') and (d') from mGW in the thermal dynamic limit at density 
no = 0.25 and no = 0.75. The left column (a), (c), (a'), (c') and the right column (b), 
(d), (b'), (d') refer to the hard-core case and the soft-core case, respectively. 
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5.2. Anyon Quasi-momentum Distribution 

According to the Eqs. (14) and (A.4) in the Appendix, the anyonic correlation function 
is written as 


^ [A -f Yl w , 

(35a) 

i<l<j 


^ [A + Y[ w*. 

(35b) 


j<l<i 


where w = = WD^. If 0 = 0, we have w = F^ + Ff + F| = 1, 

otherwise the absolute value is hh < 1. 

In the Bose limit 0 = 0, the correlation function of anyons is 

(ajaj) = dijUo + (1 - 5ij)C . (36) 


Using Eq. (15b), the ground-state quasi-momentum distribution of anyons can be 
written as 

{nt'>)=no-C + CL6k,o. (37) 

If 6* 7 ^ 0 and then |t<;| < 1, the quasi-momentum distribution of anyons is written 
as 


(hf^) = no+ 


.kiL-l)-Lu+u^ 


(1 


-c.c. 


u 


B 

Tv 


^_;k {L-l)-Lv+v^ 

{l-vY 


f c.c. 


(38) 


where u = we and v = we In the thermodynamic limit, this reduces to 


fa) \ L —>-00^ 


{nT) 


| U !|<1 


) A 

> Uq + A 


2 cos A; — 2IU cosy 
l-2IUcos(fc-x) + lU2 


+ B- 


2 cos k — 21U cos(6' — x) 
-2Wcos[{k + 0)-x] + W‘^ ’ 


(39) 


where W = \w\ and y = arg(ta). 

In the hard-core limit F^ = 0, so we have ta = z, y = yj, A = ?7,o(l — uq), and 
B = 0, thus Eq. (39) reduces to Eq. (26). For small densities no 1, which implies 
F 2 ^ ~ 0, we expect the quasi-momentum distribution of anyons in the soft-core case to 
be similar to that in the hard-core case. This can be shown by comparing Fig. 4 (a) 
and (b) from DMRG or (a') and (b') from mGW at the density uq = 0.25. However, 
for the fractional phase 0 = n in the soft-case, we also note that there exists a quite 
small area with the property > 1 that is shown as the shaded (yellow in color) 

zone in the insets of Fig. 4 (b) and (b'). The shaded area increases while increasing the 
density, for example no = 0.75 is shown in Fig. 4 (d) and (d'). This result reveals that 
anyons with the statistical phase 0 = 71 are pseudofermions, i.e. possibly more than 
one identical anyon exist in one state but they behave like a fermion while exchanging 
two anyons on different sites. Comparing the results from DMRG and mGW in Fig. 4, 
we demonstrate again that the mGW approach grasps the main features of the quasi¬ 
momentum distribution of anyons. Finally, we should note that the hard-core case is 
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- n„ =0.25 mGW 

- n„ =0.50 mGW 

- - n„ =1.25 mGW 

- - n„ =1.50 mGW 
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n„ =0.50 DMRG 
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fa) 

Figure 5. Quasi-momentum Kiaax plotted against the fractional phase 0 of anyons 
resulting from DMRG {L = 120) and mGW (in the thermodynamic limit L —>■ oo) 
with parameters J/U = 0.1, density no = 0.25,0.5,1.25, and 1.5, respectively. The 
uncertainty of the quasi-momentum is of the order Afc = 27r/T. 


more suitable to describe the step-like behaviour of anyons {6 = tt) than the soft-core 
case. 

Fig. 5 shows the fractional phase ^-dependence of anyonic iFmax at the densities 
no = 0.25, 0.5, 1.25, and 1.5. In general, for the density 0 < no < 1 and 0 < 6^ < vr, iFmax 
at density no = 0.5 is larger than that at no = 0.25, which coincides with the hnding of 
the hard-core case. Moreover, for 1 < no < 2, the DMRG results show that the peak 
position Kmlx of the quasi-momentum distribution at density no = 1.25 coincides with 
that at density no = 0.25 within uncertainty. This is also shown by the DMRG results 
at density no = 1.5 and no = 0.5. However, the mGW approach in the thermodynamic 
limit shows that iFmax at density no-l-l is a little larger than that at density no, especially 
for the zone 0 ~ tt. 

6. Conclusions 

In summary, we have studied the ground-state property of the ID Anyon-Hubbard 
model. With the help of a fractional version of the Jordan-Wigner transformation, the 
Anyon-Hubbard model is mapped to the occupation-dependent hopping Bose-Hubbard 
model and, thus, the Hilbert space of anyons can be constructed from that of bosons. 
By calculating the two-point correlation function of creation and annihilation operators 
of bosons and anyons, we investigate the quasi-momentum distributions interpolating 
between Bose-Einstein statistics and Fermi-Dirac statistics. Theoretically, in order to 
include the influence of the fractional phase of anyons on the many-body wavefunction, 
we modify the classical Gutzwiller mean-held approach and get an analytic expression for 
the quasi-momentum distributions of anyons and bosons. In order to test the accuracy 
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of the niGW mean-field approach, we use DMRG for numerical calculations. 

In the hard-core case, the results show that the bosons condensate at 
zero momentum and have a symmetric quasi-momentum distribution around zero 
momentum. Due to the nonlocal string property, more anyons are shifted to a positive 
momentum and have an asymmetric quasi-momentum distribution, where the peak 
position depends on both the fractional phase and the particle number density. 

In the soft-core case, the results show that the quasi-momentum peaks of bosons 
strongly depend on the particle number density. At density uq < 1, A/Sq = 2m7i leads to 
a quasi-momentum peak of bosons at A; ~ 0. At density 1 < Uq < 2 , A/3i = 2l7r — 9 leads 
a quasi-momentum peak of bosons at /c ~ —6. Again, the quasi-momentum distribution 
of anyons show nonlocal string behaviour and yield similar features as in the hard-core 
case. Furthermore, anyons with 6 = tt are pseudofermions, i.e., there exists more than 
one identical anyon in one state in the soft-core case. However, in the hard-core case, 
the quasi-momentum distribution of anyons is a typical step-like function. 
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Appendix A. Gutzwilier Mean-Field Approach 

Based on the Gutzwiller wave function (5), the interaction energy turns out to be 



j 


(A.l) 



j n=0 


whereas the chemical potential term within a grand-canonical description reads 



(A.2) 


In the same way, the expectation value of the kinetic energy term can be expressed as 




J 


(A.3) 




J 


(A.4a) 


where the expectation values {Oj) = yield 

^max 

(i,) = E 4''*/ihv^rn: = (i*)*, 

n =0 

^max 

, (A.4b) 

n =0 

^max 

= (A-4c) 

n={) 

^max 

(g±i0n^) _ 

n =0 

Using Eq. (A.4), the Gutzwiller kinetic energy reduces to 

^max 

4i„ = - j E E y("+1)(™+1) [tf’/S:/U''7UU"'"+c.c,]. (a,5) 

j m,,n 


Using the polar decomposition fY = Fn 


vO) 


we have 


f{j) fU)* x(i+i)* x(i+i) _ rp{j) p(j) pO'+i) pii+^)pi[YY 

Jn J n+lJ m J m+1 ^ n ^ n+l-'- m ^ m+1 i 


(A.6) 


where /9n ^ denotes the difference between the phase of the probability 

amplitude of Ending n bosons and f^h of finding n + l bosons on site j. 

The ground state of the system with given parameters is determined from varying 
the coefficients fn'^ and minimizing the total energy. This can either be done for a 
given chemical potential fi, or for a given mean particle number no = (n) = YlnYo 
according to Eq. (A.4c). In this paper, we assume that no is given, so the chemical 
potential = —pnoT resulting from (A.2) does not depend on fY\ Thus, after having 
determined the Gutzwiller coefficients fn \ the chemical potential /i is fixed by 




1 d{Sint + i^kin) 

L driQ 


(A.7) 


Appendix A. 1. Hard-Core Case 

In the hard-core limit, bosons are impenetrable and each site contains at most one 
particle, i.e., n^ax = 1- In the GW approach, the normalization condition (6) reduces 
to Fq + FI = 1. For a given particle number N, the mean density of particle number 
is given by uq = N/L = F^ according to the Eq. (A.4c). Thus the absolute values of 
both Gutzwiller amplitudes are determined according to 
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The expectation value of the energy per lattice site E^ot = ^tot/L is then given by 

- 2Jno(l - uq) , 

= -f^no - 2Jno(l - no) cos A/3o , 


(A.9) 


where E^^"^ and refer to the cGW and niGW approach, respectively. The ground 

state is determined by minimizing the energy E^^"^ and . This leads to 


A/So = 2m7r, 


(A.IO) 


where m = 0, ±1, ±2, • • •. This expression shows that the mGW approach reduces to 
the cGW approach for hard-core anyons. 

Appendix A.2. Soft-Core Case 

Using the normalization condition (6) and Eq. (A.4c) in the GW approach, we have 

F2 ^ p2 + = 1, 

El + 2E2 = no , 


(A.ll) 


where we suppose rimax = 2 as an approximation. The expectation value of the energy 
per lattice site reads now 


E^?^=-tino + UEf - 2JEf\ E^ cos(A/So) + 2Ef cos(A/Si + 9) 


n dip IP f ^Po + ^f3l + 9\ 

2v2FoF2Cos^- - - j cos 


(^Po—(3i- 


(A.12) 


In order to minimize (A.12), we must choose for the respective phases 

A/So = 2m7r, A/Si + 6 = 21 tt , 


/^o — /^i 


- = (m + OvT, 


(A.13) 


where m,l = 0, ±1, ±2, • • •. This demonstrates that the non-vanishing 6 influences the 
phase of the many-body wavefunction. Taking into account Eq. (A.ll), the ground 
state can be found by varying the Gutzwiller coefficients Eo,Ei, E 2 and minimizing the 
energy 

EZt^ = -t^no + UE^ -2JC, (A. 14) 

where C = Ef{Ef + 2y/2EoE2 + 2Ef). Note that the resulting occupations Ef, Ef, and 
Ef of the states |0), |1), and |2) are not affected by 6. Obviously, in the hard-core limit 
J/U —)■ 0, we reproduce E 2 —)■ 0. Thus, in that case Eq. (A.ll) goes over to (A.8), so 
the calculation reproduces the hard-core case of the last section. 

Furthermore, if we set A/Sq = 2 m7r and A/Si = 2171 in Eqs. (A.12), the mGW 
approach reduces to the cGW approach, for which the energy is 


77 .CGW _ 
■^tot 


-/ino-|-UF 2 ^ —2JF^^< Eq 2 E 2 cos{9) 


2 V 2 E 0 E 2 cos(^/3o “ } 


9' 


(A.15) 
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Comparing Eq. (A. 14) with (A. 15), it is straightforward to conclude that < 

which shows that the niGW approach is superior to the cGW approach. 
Therefore, we restrict ourselves to evaluate the mGW approach for the soft-core case. 
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